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SUMMARY/OVERVIEW 

This  report  describes  progress  during  the  third  year  of  a  collaborative  research  program 
combining  the  expertise  of  individuals  from  several  universities  to  develop  a  new  ability  to 
predict  the  propulsion  performance  of  solid  rocket  motors.  The  focus  of  the  research  on 
nonsteady  behavior  is  imique  and  the  overall  project  is  not  possible  at  any  one  of  the  institutions 
participating  in  this  coordinated  research.  The  individual  tasks  which  we  are  studying  will 
pursue  solid  propellant  decomposition  under  unsteady  conditions,  nonsteady  aspects  of  gas  phase 
flame  structure  measurements,  numerical  modeling  of  multidimensional  flame  structure, 
propellant/flame  interactions  and  overall  nonsteady  propellant  combustion  characteristics  in 
realistic  rocket  motor  environments.  Our  goal  is  to  develop  general  models  of  fundamental 
mechanisms  of  combustion  instability  that  can  be  applied  to  a  variety  of  new  energetic  materials. 

TECHNICAL  DISCUSSION 

The  objectives  of  MPMS  smdies  were  to  determine  the  chemical  mechanisms  for  HMX, 
GAP,  and  HMX/GAP(80%  HMX  and  20%  GAP  by  weight)  propellants  during  steady  state 
laser-assisted  burning,  and  to  determine  the  response  of  gas-phase  species  and  burning  rate  of 
HMX  during  sinusoidal  laser  heating.  A  triple  quadrupole  mass  spectrometer  (TQMS)  and 
quartz  microprobes  are  used  o  obtain  quantitative  species  profiles  above  the  surface  of  the 
sample  during  CO2  laser  induced  decomposition  or  during  laser  assisted  deflagration  of  the 
samples.  Gas  phase  temperature  profiles  are  also  measured  with  micro-thermocouples  no  larger 
than  75  micron.  The  results  of  the  studies  provide  the  boundary  conditions  at  the  interface 
between  the  condensed  and  gas  phases  which  is  critical  for  modeling.  In  addition  they  provide 
gas  phase  species  and  temperature  profiles  over  a  range  of  experimental  conditions  that  will  be 
used  to  test  the  models. 

Steady  State  Measurements 

HMX,  GAP,  and  HMX/GAP(80%  HMX  and  20%  GAP  by  weight)  propellants  were 
investigated  at  atmospheric  pressure  in  argon  with  laser  heating.  Because  most  of  the  work  has 
been  summarized  in  past  reports,  only  the  HMX/GAP  propellant  results  will  be  discussed  here. 
For  the  HMX/GAP  propellant,  a  long  ignition  delay  of  1.5  sec  was  usually  observed  at  a  heat 
flux  of  100  W/cm^,  while  an  ignition  delay  of  42  ms  was  found  at  a  heat  flux  of  300  W/cm^.  At 


a  heat  fliix  of  100  W/cm^,  a  significant  amount  of  carbonaceous  residue  was  observed  on  the 
burning  surface.  In  contrast,  a  small  amount  of  carbonaceous  residue  was  found  only  at  the  edge 
of  the  propellant  at  a  heat  flux  of  300  W/cm^. 

Figure  1  displays  species  profiles  for  HMX/GAP  at  atmospheric  pressure  in  argon  with  a 
heat  flux  of  100  W/cm^.  The  species  profiles  display  three  distinctive  regions;  the  primary 
reaction  zone,  the  preparation  zone,  and  the  secondary  reaction  zone.  The  primary  reaction  zone 
extends  from  0  to  0.5  mm  above  the  surface.  The  reactions  are  evidenced  by  the  fact  that  H2O 
and  NO2  are  consumed  as  NO,  CO,  CO2,  and  H2O  are  produced.  A  steep  temperature  rise  was 
observed  just  above  the  surface,  which  is  in  line  with  the  highly  exothermic  reaction  near  the 
surface.  At  the  end  of  the  primary  reaction  zone,  NO2  is  almost  totally  consumed  but  some 
H2CO  remains.  The  preparation  zone  lasts  from  0.5  to  2.7  mm  above  the  surface.  In  this  region, 
all  the  species  profiles  are  relatively  flat,  which  is  consistent  with  the  observed  flat  temperature 
profile  after  the  sharp  temperature  rise.  The  secondary  reaction  zone  starts  from  2.7  mm  above 
the  surface,  where  NO  and  HCN  are  consumed  and  N2  and  CO  are  produced  at  a  similar  rate. 
After  the  NO  is  completely  consumed  in  the  secondary  reaction  zone,  some  HCN  remains,  at  a 
mole  fraction  of  0.02,  which  is  consistent  with  the  fuel-rich  composition  for  the  HMX/GAP 
propellant.  It  is  noteworthy  that  the  mole  fraction  of  H2O  did  not  increase  from  3.5  to  4.2  mm 
above  the  surface.  This  is  consistent  with  the  fact  that  HCN  and  NO  will  react  without  forming 
water. 

Though  only  20%  GAP  was  added  into  pure  HMX,  the  gas-phase  zonal  structures  were 
significantly  different  between  laser-assisted  combustions  of  HMX/GAP  and  pure  HMX.  The 
HMX/GAP  flame  displayed  a  narrow  primary  reaction  zone  and  a  relatively  long  preparation 
zone  between  the  primary  and  secondary  reactions,  while  the  pure  HMX  flame  showed  a  long 
primary  reaction  zone  and  no  preparation  zone  before  the  secondary  reaction  zone.  The  high 
exothermic  condenses-phase  reaction  of  GAP  could  increase  the  condenses-phase  temperature  of 
HMX/GAP,  and  the  observed  increase  in  gas-phase  temperature  'will  accelerate  the  reaction 
between  H2CO  and  NO2  near  the  surface.  Thus,  the  primary  reaction  zone  was  considerably 
shorter  in  the  HMX/GAP  flame.  Several  factors  could  account  for  the  fact  that  the  secondary 
reaction  zone  did  not  occur  right  after  the  primary  reaction  zone  in  the  HMX/GAP  flame.  The 
higher  burning  rate  for  HMX/GAP  could  be  part  of  the  reason  because  the  higher  gas  velocity 
requires  a  longer  distance  to  achieve  the  same  reaction  time.  The  fuel-rich  gas  mixture  could 
also  be  a  contributing  factor  because  addition  of  GAP  increased  the  fuel  to  oxidizer  ratio.  It  was 
evident  that  the  mole  fraction  of  HCN  at  the  surface  for  the  HMX/GAP  flame  was  higher  than 
that  for  the  HMX  flame.  Thus,  the  higher  gas  velocity  and  the  large  amount  of  HCN  could 
significantly  delay  occurrence  of  the  secondary  reaction  in  the  HMX/GAP  flame. 

Unsteady  Measurements  for  HMX 

Over  the  last  year  of  the  program,  a  new  experimental  facility  was  designed  and 
constructed  which  permitted  the  study  of  the  response  of  HMX  to  sinusoidal  laser  heat  fluxes. 
(This  facility  was  ^ded  mainly  through  the  MURI  funds  from  ONR  and  BMDO.)  With  this 
system,  thrust,  species,  flame  height,  and  light  emission  responses  to  a  sinusoidal  heat  flux  were 
measured  at  frequencies  from  3.9  to  250  Hz  and  mean  heat  fluxes  from  35  to  90  W/cm^  for 
HMX  pellets.  The  CO2  laser  with  associated  electronics  and  software  was  used  to  provide  a 
sinusoidal  heat  flux.  The  laser  system  has  been  proven  to  be  an  effective  tool  for  the  tests  of  the 
laser-driven  response  function.  All  tests  were  done  using  a  heat  flux  sweep  with  discrete 
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Figure  1:  Species  profiles  above  the  surface  of  HMX  during  laser-assisted  combustion  at 
one  atmosphere  and  a  heat  flux  of  100  W/cm^ 

increments  of  frequency.  The  time  interval  for  each  frequency  always  covered  several  cycles. 
This  unique  sweep  technique  has  the  advantages  of  using  much  less  material,  normally  one 
single  sample,  to  cover  a  given  frequency  range  and  of  reducing  the  uncertainties  arising  from 
separate  tests.  Thus,  less  material  was  consumed,  and  relatively  good  data  was  produced. 

The  microprobe/TQMS  system  was  successfully  applied  for  measuring  the  species  under 
oscillating  conditions.  The  stable  species  response  to  sinusoidal  laser  heating  had  not  been 
obtained  by  earlier  researchers.  The  results  for  the  secondary  flame  region  of  HMX  showed  that 


NO  always  oscillated  in  phase  with  HCN,  but  180°  out  of  phase  with  N2.  The  results  were  in 
line  with  the  relations  between  product(N2)  and  reactants(HCN  and  NO). 

Figure  2  is  a  plot  of  the  HMX  oscillation  amplitudes  and  phases  for  the  flame  height, 
thrust,  light  emission,  and  species  responses  as  a  function  of  frequency.  The  four  sets  of  data 
were  obtained  at  a  mean  heat  flux  of  60  W/cm^.  The  thrust  phase  response  shows  a  phase  lead  at 
frequencies  below  12  Hz  and  a  phase  lag  approaching  100°  at  125  Hz.  Theoretically,  the 
resonant  peak  should  be  at  a  frequency  where  the  phase  changes  from  lead  to  lag.  Both  the 
thrust  amplitude  and  phase  relations  indicate  a  resonant  peak  of  1 1-1 5  Hz.  The  light  phase 
always  lagged  behind  the  thrust  phase,  and  the  thrust  and  light  phases  diverged  with  an  increase 
in  frequency.  The  species  phase  was  negative  at  all  frequencies,  but  the  amplitude  response 
suggests  at  peak  at  approximately  1  OHz.  The  difference  between  the  thrust  and  species  response 
functions  suggests  that  the  species  composition  at  the  surface  varied  with  burning  rate,  which 
was  observed  during  self-oscillating  burning,  to  be  discussed  in  the  next  paragraph. 

Oscillating  burning  was  also  observed  at  a  constant  heat  flux  of  30  W/cm^.  Simultaneous 
surface  temperature  and  species  measurements  were  performed  at  the  condition.  This  unique 
experimental  technique  enabled  the  relationship  of  the  surface  temperature  to  all  the  stable 
species  at  the  surface  under  the  transient  conditions  to  be  determined.  During  the  test, 
oscillations  of  species,  temperature,  and  regression  rates  were  observed  and  were  found  to  occur 
at  a  frequency  of  approximately  4  Hz.  NO2,  HCN,  and  triazine  oscillated  in  phase  with 
temperature,  whereas  N2O,  H2CO,  and  the  species  at  mass  28  were  180°  out  of  phase  with 
temperature.  The  production  of  NO2  and  HCN  was  favored  with  an  increase  in  temperature  and 
burning  rate,  while  the  production  of  N2O  and  H2CO  became  more  important  with  a  decrease  in 
temperature  and  burning  rate.  This  was  qualitatively  in  line  with  the  competing  reaction 
branches  proposed  by  Brill  and  others  for  the  condensed-phase  reaction  of  HMX. 

Gas  Phase  Flame  Processes  (Branch) 


The  development  of  detailed  models  of  solid  propellant  combustion  requires  detailed 
information  on  the  gas-phase  chemistry  occurring  above  the  propellant  surface.  Different  flame 
systems  have  been  identified  as  important  to  study  in  order  to  gain  greater  understanding  of  solid 
propellant  gas-phase  chemistry.  In  particular,  systems  that  need  study  are  flames  of  carbon 
monoxide  and  hydrocarbons  burning  with  nitrogen  oxides.  Since  both  premixed  and  diffusion 
flame  chemistries  exist  above  the  propellant  surface,  premixed  and  diffusion  flames  involving 
these  reactants  need  to  be  studied.  In  the  past  decade,  several  experimental  studies  have  been 
made  of  premixed  flames  consisting  of  these  reactants.  More  recently,  counterflow  diffusion 
flames  consisting  of  these  reactants  have  been  studied  as  well.  In  the  present  work,  several 
different  premixed  and  counterflow  diffusion  flames  have  been  modeled  using  a  common,  275- 
reaction  mechanism. 

Solid  propellant  combustion  models  require  gas-phase  chemical  kinetic  mechanisms  that 
are  as  small  as  possible.  Therefore,  for  each  premixed  flame  system  studied,  sensitivity  and  rate 
generation  analyses  have  been  used  to  identify  which  reactions  in  the  detailed  mechanism  are 
critical  to  each  system.  Based  on  these  analyses,  and  comparisons  with  the  experimental  data, 
reduced  mechanisms,  consisting  of  only  the  essential  chemical  reactions  for  each  flame  system, 
have  been  developed.  These  reduced  mechanisms  have  been  combined  into  a  comprehensive, 
reduced  mechanism  consisting  of  42  reactions.  In  general,  this  reduced  mechanism  models  the 
flames’  chemistry  as  accurately  as  the  full  mechanism  does. 
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Figure  2:  HMX  gas-phase  species  amplitude  and  phase  responses  compared  with  thrust, 
flame  height,  and  light  emission  signals. 


The  foil  and  reduced  mechanisms  have  also  been  used  to  model  methane-nitrous  oxide 
and  carbon  monoxide-nitrous  oxide  counterflow  diffusion  flames.  From  comparisons  with  the 
experimental  data,  it  is  determined  that,  in  general,  the  foil  mechanism  and  the  reduced 
mechanism  model  the  diffusion  flame  chemistry  as  accurately  as  they  modeled  the  premixed 
flame  chemistry.  Also,  in  the  diffusion  flames,  there  is  reasonable  agreement  between  the 
modeling  results  using  the  foil  and  reduced  mechanism. 

The  problem  of  N2O  decomposition  remains  a  key  area  requiring  further  study.  The  rate 
constants  for  this  reaction  had  to  be  modified  to  model  the  CO-N2O  flames  and  the  CH4-N2O 
diffusion  flame  accurately.  The  reason  for  this  discrepancy  is  still  to  be  determined,  but  future 
studies  of  the  N2O  decomposition  reaction  are  suggested. 

Monopropellant  Burning  Model  (Smooke) 

Our  research  has  focused  on  the  development  and  application  of  computational  models 
that  can  be  used  to  elucidate  information  related  to  the  burning  of  nitramine  propellants.  The 
work  focused  on  three  specific  projects,  each  of  which  is  described  in  some  detail  below. 

A  model  was  developed  that  applied  detailed  chemistry  and  complex  transport  to  a  two- 
dimensional  flame  in  which  a  fuel  stream  was  sandwiched  between  two  oxidizer  streams.  The 
work  focused  on  the  burning  of  foel-N20  flames  in  a  Wolfhard-Parker  configuration.  Unlike 
some  models  in  which  diffusion  in  the  downstream  direction  was  neglected,  we  have  treated  the 
folly  elliptic  problem.  A  discrete  solution  was  obtained  by  combining  a  steady-state  and  a  time- 
dependent  solution  method.  A  time-dependent  approach  was  used  to  help  obtain  a  converged 


numerical  solution  on  an  initial  coarse  grid  using  a  flame  sheet  starting  estimate.  Grid  points 
were  inserted  adaptively  and  Newton’s  method  was  used  to  complete  the  problem.  The  model 
was  used  to  compare  temperature  and  species  profiles  with  experimental  measurements  made  at 
the  University  of  Colorado.  Both  CH4-N2O  and  CO-N2O  flames  were  considered.  The  results  of 
this  investigation  enabled  the  probing  of  the  fluid  dynamic  and  thermochemistry  interaction  and 
its  effect  on  the  structure  of  fuel/N20  flames. 

A  computational  model  was  developed  which  is  applicable  to  the  study  of  HMX 
combustion.  The  three-tiered  model  allowed  for  solid,  liquid  and  gaseous  phases  of  HMX.  The 
resulting  nonlinear  two-point  boundary  value  problem  was  solved  by  Newton’s  method  with 
adaptive  gridding.  In  the  model  the  burning  rate  was  computed  as  an  eigenvalue  which  removed 
the  uncertainty  associated  with  employing  evaporation  and  condensation  rate  laws  in  its 
evaluation.  We  have  applied  the  model  to  the  study  of  both  laser-assisted  and  self-deflagration 
of  HMX  monopropellants  and  we  have  compared  the  results  with  experimental  measurements 
made  at  China  Lake.  The  burning  rates  have  been  computed  over  a  wide  range  of  pressures  and 
they  compare  will  with  experimental  results  from  one  to  ninety  atmospheres.  The  burning  rate  is 
found  to  be  proportional  to  the  pressure  raised  to  the  0.82  power.  Sensitivity  of  the  burning  rate 
to  initial  propellant  temperature  was  calculated  and  found  to  be  extremely  low  which  is  in 
agreement  with  past  theoretical  predictions  and  experimental  data.  Results  for  self-deflagration 
studies  do  not  indicate  a  two-zone  flame.  Laser-assisted  studies  illustrate  a  distinct  primary  and 
secondary  flame  separated  by  a  dark  zone,  the  length  of  which  is  dependent  upon  the  incident 
laser  flux  intensity. 

The  HMX  model  was  modified  to  study  the  combustion  of  HMX  propellants  over  the 
pressure  range  of  90-3400  atm.  New  to  this  effort  was  the  treatment  of  the  gas-phase  as  non¬ 
ideal.  The  Becker-Kistiakowsky-Wilson  (BKW  equation  of  state  was  used  to  relate  the 
compressibility  factor  to  the  pressure,  density,  gas  constant  and  temperature.  Thermodynamic 
properties  were  then  derived  from  standard-state  (ideal  gas)  mixture  properties  plus  real  gas 
mixture  departure  functions.  Reverse  rate  constants  were  obtained  from  the  forward  rate 
constants,  compressibility  factor  and  the  fugacity  based  equilibrium  constants.  Results  of  the 
deflagration  rate  as  a  function  of  pressure  were  found  to  be  in  good  agreement  with  recent 
experimental  data.  Differences  in  HMX  bum  rates  predicted  by  the  ideal  and  non-ideal  gas- 
phase  models  are  attributed  to  the  resulting  heat  flux  differences  at  the  HMX  gas-liquid  interface. 

Propellant  Combustion  and  Chamber  Dynamics  lYang) 

Results  obtained  from  the  preceding  tasks  provide  a  concrete  basis  for  understanding  the 
detailed  combustion  mechanisms  of  solid  propellants  in  an  isolated  environment.  To  extend 
these  fundamental  investigations  to  practical  motor  design  issues,  it  is  important  to  study  the 
propellant  burning  behavior  under  conditions  representative  of  realistic  rocket  motor 
environments.  Most  previous  research  on  motor  interior  ballistics  was  either  focused  on  cold- 
flow  simulation  or  based  on  simplified  approaches  with  propellant  burning  rates  modeled  by 
empirical  formulas.  Very  limited  is  known  about  detailed  physico-chemical  processes  near  the 
propellant  surface  which  dictate  the  burning  characteristics.  The  purpose  of  this  task  is  to 
remedy  these  technical  and  scientific  deficiencies  by  treating  the  intricate  interactions  between 
local  flow  motions  and  chemical  reactions  at  scales  sufficient  for  resolving  the  rapid  variations  of 
flow  properties  in  the  near  field.  Both  steady  and  oscillatory  conditions  are  considered.  The 


work  mainly  addresses  two  issues:  how  the  motor  internal  flow  is  established  by  propellant 
burning,  and  how  the  local  flowfield  affects  propellant  burning. 

During  the  past  year,  a  series  of  analyses  have  been  conducted  to  study  the  combustion 
characteristics  of  AP/HTPB  composite  propellants  in  two-dimensional  chambers,  as  shown  in 
Fig.  3.  The  formulation  is  based  on  the  complete  conservation  equations  for  both  gas  and 
condenses  phases,  with  full  account  taken  of  variable  transport  and  thermodynamic  properties. 
Only  the  reduced  chemical  kinetics  schemes  established  in  the  flame  chemistry  studies  are 
incorporated  into  the  flow  solver  to  render  numerical  calculations  manageable.  Turbulence 
closure  is  achieved  using  a  modified  turbulent  transport  model  which  takes  into  account  the  wall¬ 
blowing  effect  arising  from  propellant  burning.  Specific  processes  investigated  include: 

•  Motor  internal  flow  development  over  a  broad  range  of  operating  conditions,  including  onset 
of  turbulence  and  its  ensuing  growth; 

•  Detailed  flame  structure  near  the  propellant  surface,  including  the  micro-scale  energy 
cascade  from  exothermic  reactions  to  flow  motions; 

•  Coupling  between  gas  and  condenses  phases,  including  subsurface  degradation  and  surface 
regressions  processes; 

•  interactions  between  acoustic  wave  and  flame  dynamics  in  the  gas  phase,  including  the 
mutual  couplings  among  the  acoustic,  vortical,  and  entropy  modes  as  well  as  their  influences 
on  heat-release  mechanisms,  and 

•  dynamics  responses  of  propellant  burning  to  local  flow  oscillations,  including  transient 
variations  of  subsurface  heat  release  and  thermal-wave  penetration. 

The  model  appears  to  be  the  first  and  most  complete  of  its  kind  to  date,  providing  comprehensive 
information  about  propellant  combustion  and  motor  chamber  dynamics.  The  instantaneous 
burning  rate  of  propellant  is  calculated  as  part  of  the  solution  from  first  principles.  Figure  4 
shows  a  typical  result  of  the  spatial  distribution  of  Rayleigh’s  parameter.  The  problem 
investigated  involves  a  two-dimensional  chamber  loaded  with  composite  propellant  subject  to 
periodically  forced  oscillations  at  the  chamber  exit.  This  result  can  be  effectively  used  to 
identify  the  day  mechanism  responsible  for  driving  combustion  instabilities  in  rocket  motors.  In 
particular,  the  effects  of  spatial  distribution  and  tempor^ll  evolution  of  the  heat-release  process,  as 
well  as  their  mutual  coupling  with  local  flow  motions,  on  motor  stability  can  be  clearly 
examined. 


Unsteady  Motor  Flows  Coupled  with 


Fig.  3  Schematic  diagram  of  a  rocket  motor  with  AP/HTPB  composite  propellaht. 


Fig.  4.  Spatial  Distribution  of  Rayleigh’s  Parameter,  <pQ*: 
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